Constraining freeze-out with yields and fluctuations 
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Abstract. We show that the simultaneous measurement of yields and fluctuations in heavy 
ion collisions is capable of falsifying and constraining the statistical hadronization model. We 
^— ^ , show how such a measurement can test for chemical non-equilibrium, and distinguish between 

■ a high temperature chemically equilibrated freeze-out from a supercooled freeze-out with an 

over-saturated phase space. We further explain how this measurement can be used to obtain a 
, model-independent estimate of the difference between "chemical" and "thermal" freeze-out. 

o 

2h : 

The statistical hadronization model (SHM) p^ElE] has been extensively applied to the study 
of soft particle production in hadronic systems. When it includes the full resonance spectrum 4 , 
the SHM can describe quantitatively the yields of all hadrons produced in heavy ion collisions. 

The ability of the SHM to describe not just averages, but event-by-event multiplicity 
fluctuations has however not been widely investigated, and its applicability is currently a matter 
^ ■ of controversy. Event-by-event fluctuations are subject to intense current theoretical and 

experimental interest 013) as a constraint for existing models and as a signature of new physics. 

This study illustrates the use of both yields and fluctuations as a probe capable of constraining 
the SHM and differentiating between freeze-out scenarios. 

The statistical hadronization model assumes that particles are produced according to a 
probability determined by their phase space density. The first and second cumulants of this 
probability distribution give, respectively, the average value (over all events) of the desired 
observable, and its event-by-event fluctuation. 

Conserved quantities can be treated in several ways, appropriate to different experimental 
situations: If the totality of the system is observed, than conserved quantities can not fluctuate. 
If a small fraction equilibrated with the rest of the system is observed, than conserved quantities 
will fluctuate event-by-event. Rigorous conservation is known as the Microcanonical ensemble, 
while allowing energy and other conserved quantities to fluctuate between the system and the 
bath leads, respectively, to the Canonical and Grand Canonical (GC) ensembles. All fluctuations 
are ensemble-specific even in the thermodynamic limit |S]. 

In this work, we use the GC ensemble, implemented in open-source software jH] to calculate 
fluctuations and yields. We motivate this choice by the fact that RHIC experiments observe the 



mid-rapidity slice of the system, comprising roughly 1/8 of the total multiplicity, an appropriate 
fraction for a GC prescription. Boost invariance, a good symmetry around mid-rapidity, links 
this rapidity slice with a sector in configuration space. If the system observed at RHIC is a nearly 
ideal fluid, the matter created in this space should be in equilibrium, grand-canonically, with 
the unobserved regions. If freeze-out temperature throughout observed space is approximately 
constant, the GC ensemble should be able to describe both yields and fluctuations |10[ 111). 

The final state yield of particle can then be computed as a function of the particle mass and 
resonance decay tree, as well as temperature and chemical potentials. (The technical details of 
the calculation are given in a parallel paper |12|). The chemical potential, implemented here via 
the fugacity A = , is the main distinguishing feature between competing freeze-out models. 
Provided the law of mass action holds, it should be given by the product of charge fugacities 
(flavor, isospin etc.). It is then convenient to parametrize it in terms of equilibrium fugacities 
A eq and phase space occupancies 7. For a hadron i with q(q) light quarks, s(s) strange quarks 
and isospin ^3 the fugacity is then 

Ai = ^+H + ~ S , Af = A|-?Ar^ (1) 

If the system is in chemical equilibrium then detailed balance requires that j q = j s = 1. In an 
expanding system, however, the condition of chemical equilibrium no longer holds. Kinetically, 
this occurs because collective expansion and cooling will make it impossible for endothermic and 
exothermic reactions, or for creation and destruction reactions of a rare particle, to be balanced. 
Provided the system remains in local thermal equilibrium, Aj can still be used as a Lagrange 
multiplier for the particle density, and the first and second cumulants can be calculated from 
the partition function the usual way ^2lE]- However, in this case in general j q ^ 1,7 S 7^ 1. 

This picture becomes particularly appropriate if the expanding system undergoes a fast 
phase transition from a QGP to a hadron gas (HG). In this case, chemical non-equilibrium JH] 
and super-cooling ^S] can arise due to entropy conservation: By dropping the hadronization 
temperature to ~ 140MeV and oversaturating the phase space above equilibrium (7^ ~ 1.5, 7 S ~ 
2), it is possible to match the entropy of a HG with that of a QGP |15j . 

Fits to experimental data at both SPS and RHIC energies indeed support these values of 
7 9jS when these parameters are fitted. Moreover, best fit 7 9)S > 1 arises for a critical energy ^7] 
(corresponding to the energy of the K/ir "horn" and system size ^M,, as expected from the 
interpretation of j q as a manifestation of a phase transition. However, the fits performed in |17j 
have not been able to rule out equilibrium models (at SPS and RHIC the difference in statistical 
significance between equilibrium and non-equilibrium is ~ 20%), which are usually preferred for 
their smaller number of parameters to fit |2U1 121j . Equilibrium freeze-out temperature varies 
between fits, ranging from 155 to 177 |21| MeV. 

Both scenarios are physically reasonable, can describe the data, and would be instrumental in 
our understanding of hadronic matter if proven correct. In particular, the HBT puzzle suggests 
we lack understanding of the last stages of the fireball evolution. Non-equilibrium is useful in 
this respect, since it affects both system volume |221 117j and emission time |1H| . 

The reason both equilibrium and non-equilibrium are compatible with data is that in a 
fit to yields the non-equilibrium phase space occupancies j s and 7^ correlate with freeze- 
out temperature JJj) making a distinction between a T = 170 MeV equilibrated freeze-out 
(7 g = 1,7s < 1) scenario and a supercooled scenario where 7q jS > 1 problematic. 

A related ambiguity is the difference between chemical freeze-out (where particle abundances 
are fixed) and thermal freeze-out (where particles cease to interact). Equilibrium models 
generally assume a long phase between these two points, which would alter considerably the 
multiplicity of directly detectable resonances. In a Non-equilibrium supercooled freeze-out, on 
the other hand, it is natural to assume that particle interaction after emission is negligible [.15 1. 



Once again, a reliable way to probe the extent of the reinteraction would be instrumental for 
our understanding of how the fireball produced in heavy ion collisions breaks up. 

We have recently shown that event-by-event fluctuations can be used to solve the 
dilemmas discussed above. The equations in ^2] make it clear that the dependence of 
fluctuations on T and 7 9 is different, allowing us to decouple these two variables, a higher 
temperature tends to decrease fluctuations w.r.t. the Poisson value expected from Boltzmann 
statistics , since it introduces greater correlations due to an increased resonance contribution. 
Increasing 7^ will rapidly increase fluctuations of quantities related to pions, due to the fact that 
at j q > 1 rapidly approaches e miT ^ T , giving fluctuations an extra boost w.r.t. yields [TT1 IT2| . 
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Figure 1. (color online) Comparing a fluctuation and a particle yield sensitive to T and j q . 
Dotted black lines refer equal temperature (T=100-200 MeV), while long-dashed red lines refer 
to 7q (7^ = 0.7 — 1.8). Left: vs. a 7T -/ 7T +. Center: H/^vs. v(Q). Right: p°/ir~ vs. cr 7r -/ 7r +, 
which also probes the extent of hadronic interactions after chemical freeze-out. As shown in 
|llj . parameters other than T and j q do not impact the observed quantities significantly 

In we used a measurement of the charge fluctuation measure v(Q) = (AQ 2 ) / (iV c h) where 
(AQ 2 ) is the fluctuation in net charge and is the charged particle multiplicity. For the ratio, 
we used A/K~, corrected for H and cj) feed-down. 

Fig. ^ applies this procedure again, this time with the easier to measure S/</> ratio, which also 
depends mainly on T and 7^ at RHIC (and on A s at SPS). The central panel shows E/(p plotted 
against v(Q),and the left panel shows E/0 plotted against the event-by-event fluctuation in the 
ir + /ir~ ratio (labeled as o~ w +/ w -). As the left panel shows, 7 9 correlates the ir + /ir~ fluctuation 
with leaving this combination of observables dependant on the freeze-out temperature only, 
and constraining the region allowed in the SHM parameter space to a narrow band. 

It should be underlined that diagrams such as those in Fig. ^ahow f° r a test of the physical 
validity of the SHM, since the SHM requires that all yields and fluctuations be described by the 
same "f q and T. If the measurement corresponding to the left panel of Fig ^ is made, and the 
result is not in the narrow band given in the figure, or if the temperature obtained in two of the 
three panels of Fig. ^ is not the same, we can conclude that physics beyond the SHM plays a 
role. Furthermore, the nature of the deviation gives a hint to its physical origin. 

In particular, comparing fluctuations to directly detected resonances probes the interval 
between chemical and thermal freeze-out. Consider, for example, the tt + /tt~ fluctuation. The 
top and the bottom terms in this ratio are linked by a large correlation term due to the p° 
decay. This correlation probes the p° abundance at chemical freeze-out, since subsequent 
rescattering/regeneration does not alter the fact that the p° decay produced a 7r + and a ir~ . On 
the other hand, a direct measurement of the p° /ir~ ratio through invariant mass reconstruction 
measures the p° abundance at thermal freeze-out, after all rescattering. Hence, comparing the 



7r + /7r~ fluctuation to the p°/vr _ ratio provides a gauge for effect of the hadronic reinteraction 
period on particle abundances. 

The application of this method is outlined in the right panel of Fig. ^ where the measured p° 
abundance is compared to the 7r + /ir~ fluctuation. Again, the assumption of freeze-out with no 
reinteraction correlates the two observables to a narrow band Dependant on only the freeze-out 
temperature. While a model-independent quantitative prediction of eventual deviations from 
the narrow band is difficult to obtain a priori, one can infer qualitatively the likely origin of 
such deviations: A shift down the p°/vr _ axis would signal a long re-interacting phase which re- 
equilibrates the directly detectable p° to a lower "thermal freeze-out" temperature, suppressing 
their detectable yield but maintaining the correlation of the decay products. A fall in yields 
together with a rise in fluctuations would be evidence of p° melting in-medium, while a rise 
in yields together with a fall in fluctuations would mean the p° abundance is augmented, 
presumably by an in-medium mass decrease. A rise in both p° yield and fluctuation would 
be very problematic to explain in a model where statistical mechanics plays a role. 

In conclusion, we have shown that the two statistical hadronization scenarios shown to be 
applicable at SPS and RHIC, one with a chemical freeze-out at 170 MeV and a long reinteraction 
phase, the other with an explosive non-equilibrium transition from a high-entropy phase at 140 
MeV, give definite and very different predictions for the interdependence of particle yields and 
event-by-event fluctuations, allowing us to falsify either of these scenarios when both yields and 
fluctuations are taken into account. We eagerly await more published data in both yields and 
fluctuations to determine weather the non-equilibrium model is really capable of accounting for 
both yields and fluctuations in all light and strange hadrons produced in heavy ion collisions. 
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